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ABSTRACT 

We have identified four new black liole candidates in M31 globular clusters using 123 Chandra, 
and 4 XMM-Newton observations of the M31 central region. The X-ray source associated with 
Bo 163 (XB163) is a recurrent transient with the highest observed luminosity ^1.4x10'^* erg s~^, 
considerably brighter than any outbursts from neutron star transients Aql X-1 or 4U 1608-452; the 
outburst apparently started ~45 days earlier than the observed peak, hence the luminosity may 
have been considerably higher. We identified XB082, XB153 and XB185 as BHCs by observing 
"low state" emission spectra at luminosities that exceed the threshold for neutron star binaries. 
The probability that these are neutron star systems with anisotropic emission beamed towards us 
is < 4 X 10^**, and their variability suggests emission from a single source. We therefore conclude 
that these systems likely contain black holes rather than neutron stars. We have now identified 
4 persistently bright BHCs in the region; the probability that these are all background AGN is 
<lxlO^^". According to theory, the donors could be tidally-captured main sequence stars, or 
white dwarves in ultracompact binaries. We find that GCs that are particularly massive (XB082) 
or metal rich (XB144) can host bright X-ray sources in addition to GCs that are both (XB163). 
Our method may reveal BHCs in other bright X-ray sources. 



Subject headings: x-rays: general 
individual — black hole physics 

Introduction 



x-rays: binaries — globular clusters: general — globular clusters: 



The first identification of a black hole X-ray 
bin ary in a globula r clus ter (GC) was made 
by Im accarone et al.l (|2007l ). when they discov- 
ered a variable GC X-ray source in the giant 
elliptical galaxy NGC4472 with an X-ray lumi- 
nosity '^4x10^^ erg s^^; its strong variability 
ruled out blended emission from multiple X- 
ray sources. A further 3 GC black hole can- 
didates (BHCs) have since been identified by 
their high lumino sities and variability: a second 
one in NGC4472 ("Maccarone et al.l l20lil ). one in 
NGC3379 (JBrassington et al. 2010:), and one in 
NGC1399 (jShih et al.ll2010l) . A further uhralumi- 



Irwin et al. 



nous X-ray source in a NGC1399 globular cluster 
and [ N H] emission lines, prompt- 
(12010 ) to suggest that it is an inter- 
mediate mass black hole feeding from a disrupted 
star. However, this scenario could not explain 
why the nitrogen line was stron g er tha n the oxy- 



exhibits fO HI 
ing 



gen line; iMaccarone fc Warner! (|201ll ) proposed 



an alternative scenario where the X-ray source il- 
luminates the wind from a R Coronae Borealis- 
like (RGB) star; in this scenario, the black hole 
need not be of intermediate mass, and this system 
may be another black hole binary. 



1 



1.1. Black hole evolution in globular clus- 
ters 

The prior lack of observed black hole bina- 
ries in GCs greatly interested the modeling com- 
munity. It has long been thought that each 
globular cluster will form significant numbers 
of black holes; however, they are likely to first 
become concentrated in the central region of 
the cluster, then ejected through super-elastic 
enc ounters with r nulti-black hole systerns (see 
e.g. Spitzeil 1969 : Sieurdsson fc Hernquist 1993 : 
Portegies Zwart fc McMillanI l2000h " The black 
holes will be considerably more massive than the 
remaining main sequence stars; if the ratio of to- 
tal black hole mass to total main sequence mass 
exceeds a critical value, then equipartition is im- 
possible, and the black holes contract into the 
center of the cluster, interact and are expelled 
()Spitzeijll969h . 



However, iMackev et al.l (|2008l ) have found evi- 
dence for black hole retention, in the variation in 
core radius with age; the cores of young globular 
clusters exhibit a strong correlation between ra- 
dius and age, while the older clusters show great 
variation in core radius. This behaviour is seen 
in such diverse environments as the Magellanic 
Clouds, the Milky Way, and the Fornax and Sag, - 
itarius dwarf galaxies; hence iMackev et al.l ( 20081 ) 
concluded that some internal mechanism was re- 
sponsible. They ran N-body simulations with 10^ 
stars, and found two mechanisms for inflating the 
core radius. Mass loss from the massive stars 
drives the core expansion for the first few hundred 
million years, until the last supernovae. However, 
the simulations suggest that the later core expan- 
sion occurs after the first binary black holes are 
formed in the cluster center, when black hole scat- 
tering and ejection heats the core; this expansion 
lasted for the remai ning time in the simulations. 
Mackev et al. (l2008l) concluded that some black 
holes must be retained by the cluster; as black 
holes are ejected, the black hole density decreases, 
and hence scattering and ejection events become 
increasingly rare. The observed variation in core 
radii for old globular clusters may therefore be ac- 
counted for by different levels of black hole reten- 
tion. 



lar clusters. They found that black hole binaries 
formed through exchange would result in widely- 
separated, transient sources; they estimated the 
duty cycle to be <0.001, consistent with the ob- 
servations at that time. They also expected that 
black hole binaries formed by tidal capture of a 
main sequence star would be tightly bound and 
persistently bright, and inferred from the absence 
of such sources that the secondary was tidally dis- 
rupted during the capture. However, they noted 
that none of their findings ruled out black hole bi- 
naries in any GC, indeed every GC could contain 
~ 1 stel lar mass black hole. More r ecent simula- 
tions bv lMoodv fc Sieurdsson suggest that 
high metalicity clusters are more likely to retain 
black hole binaries, keeping ~5%, compared with 
1% for low metallicity clusters; the ejection mech- 
anism was more efficient for binaries with mass ra- 
tios >0.3, and higher metalicity clusters produced 
more binaries with low mass ratios. 
While iKalogera et al. 



()2004( ) only considered 



non-degenerate donor stars, GC black hole bi- 
naries could be ultra-compact systems accreting 
from white dwarfs; the formation and evolution 



of suc h systems are discussed by llvanova et al 
(|2010l ). They consider tidal capture of the WD 
to be too inefficient, and also discount the pos- 
sibility that primordial binaries will survive; in- 
stead, they consider exchange, collisions with red 
giants (the method for forming ultra-compact neu- 
tron star -|- white dwarf binaries), and three-body 
formation (which is expected to result in merg- 
ers for non-degenerate stars). Their conservative 
formation scenarios could account for the recent 
GC BHCs if ~10% of all formed black holes are 
retained, and interact with other stars as well as 
black holes; only 1 % retention is needed in their 
optimistic scenarios. iMaccarone et al. I ()2010l) sug- 



gested that the first GC black hole candidate is a 
black hole -I- white dwarf hierarchical triple, as it 
would explain the strong [ O HI] emission an d lack 
of B aimer lines found by (jZepf et al. 2008j), and 
also the strong X-ray variability. 

1.2. Properties of globular clusters host- 
ing strong X-ray sources 



Kalogera et al.l (|2004l ) explored the possibilities 



of forming black hole X-ray binaries in globu- 



Silk fc AronsI (jl975[ ) examined four X-ray bright 
Galactic GCs, and noted that their central escape 
velocities (~30-70 km s~^) were larger than for 
the majority of GCs; 30% of the GCs had escape 
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velocities >30 km s~^, while only 15% had escape 
velocities >40 km s""'^. Hence, those X-ray bright 
GCs were more massive than most. Furthermore, 
they discovered a correlation between the X-ray 
luminosities and metalicities of these clusters, with 
the least luminous X-ray source located in a metal 
poor GC, and the brighter X-ray sources in metal 
rich GCs, based on their integrated spectral types. 
Out of 83 globular clusters, only 6 had escape ve- 
locities >40 km s^^ and spectral types later than 
GO; two of these GCs contained the most luminous 



ray sources ( 


Silk & Arond 19751) 


Bellazzini et al. ( 


1995) later ( 



confirmed that 

GCs hosting X-ray sources >10'^^ erg s~^ were sig- 
nificantly denser. They found that Kolmogorov- 
Smirnov testing gave a 99.7% probability that 
the X-ray bright Galactic GCs were drawn from 
a different density distribution than those GCs 
without X-rays; similarly, the probability that 
the M31 X-ray GCs were significantly denser that 
the non-X-ray GCs was 98.8%. Since the X-ray 
bright Galactic GCs are mostly in the disc, it 
was necessary to demonstrate that the metalicity 
was a driver of the X-ray lumino sity, rather than 
a byp roduct of their locations. iBellazzini et al.l 
(|1995[ ) found a similar metalicity enhancement for 
X-ray bright GCs in M31, over non- X-ray GCs 
in the same region of the sky, hence they showed 
that metalicity, rathe r than location, was impor- 
tant. Recently, Peac ock et al.l ( 2010bl ) conducted 
a much-expanded survey of M31 GCs, associat- 
ing X-ray sources with 41 out of 416 old clus- 
ters. Their K-S tests showed a 98% probability 
that GCs hosting XBs are redder than the gen- 
eral population, and a 91% probability that they 
are more meta l rich th an the general population. 
Peacock et al.l (|2010bl ) note that the metallicity 
uncertainties are large, and likely weaken any gen- 
uine relationships between the clusters. 



When iKundu et al.l (|2002f) surveyed the X-ray 
sources in the giant elliptical galaxy NGC 4472, 
they found that ^40% of the bright X-ray sources 
were associated with GCs, a significantly higher 
fraction than for the Milky Way. However, the 
fraction of X-ray bright GCs was ^4%, remark- 
ably similar to the fractions observed in the Milky 
Way, and other galaxies studied; they therefore 
concluded that the formation of X-ray binaries in 
GCs must be infiuenced more by the GC proper- 
ties than the galaxy. They also used the strongly 



bimodal metalicity distribution to show that X-ray 
binaries were 3 times more numerous in the metal- 
rich GC population than the metal-poor GCs. 

1.3. Our method for identifying black hole 
XBs in globular clusters 

The most secure method of identifying black 
hole binaries requires the mass function to be de- 
rived from the radial velocity curve of the com- 
panion star. However, this method is unfeasible 
for X-ray binaries in globular clusters; these sys- 
tems are most likely to be located in the cluster 
center, making identifying the counterpart impos- 
sible with current instrumentation. Hence alter- 
native indicators are needed to identify GC X-ray 
sources containing black holes. 

In summary, our method for identifying black 
holes in extragalactic X-ray binaries involves 
searching for high luminosity Comptonized emis- 
sion, common to neutron star and black hole bi - 
naries at low accretion rates ( van der Klid 1994 ) ; 
any thermal component is too faint to observe 
in extragalactic X-ray binaries. At higher ac- 
cretion rates, black hole binaries can exhibit a 
high/soft state (>90% thermal emission), or a 
two-component steep power law state consisting 
of t hermal and Comptonized com ponents (see 
e.g. iMcClintock fc RemillardI l2006l . and refer- 
ences within). The various neutron star binary 
behaviors at higher luminosities all include a 
stro ng thermal contribu t ion to the emission (see 



White et al. 



19891 : Ivan der Klis 



9861: iHasinger fc van der Khi 
19941 ). A comprehensive sur- 



vey of Galactic neutr on star X-ray binaries by 
Gladstone et al.l ( 2007 ) revealed that no neutron 
star X-ray binary exhibited low state behaviour 
at 0.01-1000 keV luminosities >10% of the Ed- 
dington limit (£Edd); for a 2 Mq neutron star, 
0.1 LEdd ~2.6xl0'^'' erg s~^. Hence, if we see 
Comptonized emission (represented by a power 
law with slope ~1.4-2) and no significant thermal 
emission at luminosities significantly higher than 
3x10^^ erg s~^, then the accretor is a candidate 
black hole. We note t hat our criter i on is s imilar 
in principle to that of Barret et al. ( 2000[ ). who 
found that only black hole binaries can produce 
20-200 keV luminosities exceeding '^1.5x10'^^ erg 



Barnard et al.l ( 20081 ) found the first such black 



hole candidate in the M31 GC Bo 45, named 
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following the Rev i sed Bologna Cat alogue V4 
(|Galleti et all I2004L 120061 120071 l2009l hereafter 
RBC). The X-ray source associated with Bo 45 
exhibited low state behaviour at a 0.3-10 keV lu- 
minosity of 2 X 10'^* erg s~^, 6 times brighter 
than the low state luminosity threshold for neu- 
tron stars; it appears to have been persistently 
bright for the last ~30 year s, and is consistent 



with the predictions of .Kaloeera et al.l (120041) for 
a bina ry formed by tidal capture. iBarnard fc Kolbl 
()2009l) bund a second black hole candidate asso- 
ciated with the M31 GC Bo 144; this also appears 
to be persistently bright. 

It has long been known that the luminosities of 
state transitions vary considerably between X-ray 
binaries, and even between outburst of a single 
system; furthermore, transients exhibit hystere- 
sis where the transition from high to low state 
occurs at a lower luminos ity than the transition 
from low to hi g h state (i Mivamo to et ah 19951 
Maccaronell2003l ). iTang et al. (2010i ) recentlv con- 



ducted a survey of state transitions in Galactic X- 
ray binaries by comparing their fluxes in the the 
RXTE/ASM and Swift/BAT; a BAT to ASM flux 
ratio >1.0 indicated a low/hard state, while a ra- 
tio <0.2 indicated a high/soft state. They iden- 
tified 128 hard to soft transitions in 28 systems, 
consisting of 20 neutron star low mass X-ray bi- 
naries (LMXBs), 7 black hole LMXBs and 1 high 
mass X-ray binary. They found that the transi- 
tion luminosity was strongly correlated with the 
peak luminosity of the soft state, suggesting that 
the luminosity of transition from hard to soft state 
is governed by non-stationary parameters, such as 
the accretion rate history, rather than the accre- 
tion rate itself; however all transitions occurred 
at 15-50 keV luminosities -0.001-0.1 ^Edd (1-0 
erg s~^ in the 15-50 keV band = 1.3 erg s~^ in 
the 0.3-10 keV band assuming power law emission 
with photon index 1.7). 

We have conducted a variability study for 40 
X-ray sources associated with objects in the RBC, 
using 123 Chandra observations (Barnard et al., 
in prep); these include 31 confirmed globular clus- 
ters (GCs), 4 candidate GCs, 2 stars, 2 galax- 
ies and a HII region. We produced long term 
lightcurves from these data, and used them to 
identify possible BHCs (BHCs), looking for hard, 
non-thermal spectra at high luminosities , as i n 
iBarnard et all (|2008l ) and lBarnard fc Kolbl (|2009l ) . 



We used archival XMM-Newton observations to 
obtain high quality spectra of these sources where 
possible. We have identified 4 new BHCs as- 
sociated with the M31 globular clusters Bo 82 
(XB082), Bo 153 (XB153), Bo 163 (XB163) and 
Bo 185 (XB185). 

The X-ray source associated with the M31 GC 
Bo 163 (hereafter refe rred to as XB163) is a known 



recurrent transient. iTrudolvubov fc Priedhorskv 



( 20041 ) found it bright in one Chandra, and sev- 
eral ROSAT observations, but found no detection 
in any other Chandra observation, nor in any of 
the four XMM-Newton observations available at 
the time. Assuming an absorbed power law emis- 
sion model, with photon index 1.7 and absorption 
equivalent to 7x10^° H atom cm^'^, they obtained 
luminosities —10^^-10^® erg s~^. 

In Section [2] we discuss the observations and 
data analysis; next, we present our results in Sec- 
tion |31 including a long-term lightcurves and anal- 
ysis of the best available spectra; finally, we dis- 
cuss our findings in Section HI 

2. Observations and data analysis 

The central region of M31 has been observed 
with Chandra on a ~ monthly basis for the last 
—11 years in order to monitor transients. We have 
analyzed 78 ACIS observations and 45 HRC ob- 
servations, in order to discern the variability of 
X-ray sources associated with GCs in this region. 
We determined the position of each source from a 
merged ACIS image, using the IRAF tool imcen- 
TROID; locations of the 4 targets are provided in 
our full survey (Barnard et al., in prep). 

For each Chandra observation, we obtained 
source and background lightcurves and spectra in 
the 0.3-7.0 keV band from circular regions; the 
background region was the same size as the source 
region, and at a similar off-axis angle. The PSF 
at each source was somewhat spread, due to their 
off-axis angles; we used extraction radii of 10" for 
XB153 and XB185, 15" for XB163 and 20" for 
XB082. 

For ACIS observations, response matrix and 
ancillary response files were made. We initially 
estimated the conversion from flux to luminos- 
ity by assuming a power law emission spectrum 
with photon index 1.7, with Nu — 6.7x10^° atom 
cm""'^, then determining the unabsorbed 0.3-10 
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keV luminosity equivalent to 1 count s~^. Af- 
ter correcting for the exposure, vignetting and 
background, multiplying the source intensity by 
this conversion factor gave the source luminosity. 
Source spectra with >200 net counts were freely 
fitted. 

For HRC observations, we used only events with 
PI 48-293 to reduce the instrumental background. 
We used the WebPIMMS tool to find the unab- 
sorbed luminosity equivalent to I count s^^, as- 
suming the same emission model as for the faint 
ACIS observations. We created a 1 keV exposure 
map for each observation, and compared the ex- 
posure within the source region with that of an 
identical on-axis region, in order to estimate the 
necessary exposure correction. We multiplied the 
background subtracted, corrected source intensity 
by the conversion to get the 0.3-10 keV luminosity. 

We created long term 0.3-10 keV lightcurves for 
each source, using the luminosities obtained from 
each observation as described above. We only in- 
cluded observations with net source counts > 
after background subtraction. We fitted each long 
term lightcurve with a line of constant intensity, 
in order to ascertain the source variability. 

All four targets were observed with XMM- 
Newton; XB163 was caught in outburst once in 
2006, December, and twice in 2008, February. We 
obtained source and background spectra in the 
0.3-10 keV band from the pn instrument; the 
background regions were circular, with the same 
size as the source region, on the same CCD, and 
at similar offset angles. The appropriate response 
matrix and ancillary response file was obtained for 
each source. We neglected the MOS detectors to 
avoid pileup. 

3. Results 

We present a merged, 0.3-7 keV ACIS image of 
the central region of M31 in Fig [TJ our BHCs are 
circled and labeled. We note that XB163 looks 
rather faint in this merged image, as it is a re- 
current transient, and quiescent for much of the 
time. These merged observations cover an approx- 
imately circular region with 20' radius. We de- 
tected 433 X-ray sources in this region, and found 
428 globular clusters from the RBC. Therefore we 
expect chance coincidences of X-ray sources within 
1" of the GC centers for 0.13 out of 428 GCs. 



Figure [2] shows the 0.3-10 keV long term 
lightcurves of our four BHCs; ACIS and HRC 
observations are represented by crosses and cir- 
cles respectively. Superposed on each lightcurve 
is the best fit line of constant intensity; we give 
the luminosity, uncertainties and I degrees of 
freedom (dof) in the caption. 

Spectra were grouped to have a minimum of 20 
counts per bin, and analyzed with XSPEC version 
12.6. Each spectrum was fitted with four models 
that were representative of various spectral states 
exhibited by LMXBs (see e.g . Ivan der Klis.lQOl: 
McClintock fc Remillardll2006[ l . The low accretion 
rate state common to neutron star and black hole 
LMXBs is dominated by inverse Compton scatter- 
ing of cool photons on hot electrons, and low state 
XBs have electron temperatures ^100-300 keV 
(see e.g. iMcClintock fc Remillard .2006.'): it is well 
represented by power law emission in the 0.3-10 
keV band, with F ~1.4-2. We fitted each spec- 
trum with two models to represent Comptoniza- 
tion: a simple power law model, and a COMPTT 
model. The black hole high state was represented 
by a disk blackbody ( McClintock fc RemillardI 
20061) . A blackbody -f- power law model is rep- 
resentative of high accretion rate neutron stars; 
we expect a blackbody with kT ~l -2 keV, con- 
tribu ting ^^10-50% of the flux (see e.g. lWhite et al 



19861 ). Each model is attenuated by line of sight 
absorption. 

The properties of each spectrum are provided 
in Table [1] assuming a power law emission model. 
For each spectrum we indicate the observation 
used, and the number of net source counts; we 
then give the column density, power law index, 
X^/dof and 0.3-10 keV luminosity for the best fit 
power law model; uncertainties are quoted at the 
90% confidence level. 

3.1. XB082 

XB082 was observed in 13 ACIS and 35 HRC 
observations, due to its high off-axis angle, but 
appears to be persistently bright. Its 0.3-10 keV 
luminosity varied by a factor ^^3, up to ^3.2 x 10'^® 
erg s~^. It was also observed in the 2001 XMM- 
Newton observation of the M31 core, with a good 
exposure time of ^27 ks. 

The power law emission model yielded the best 
fit to the XMM-Newton pn spectrum; A^h (absorp- 
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tion) = 3.9±0.5 xlO^-^ atom cm ^, substantially 
higher than the Galactic absorption along the line- 
of-sight (iVn = 6.7x102" atom cm'^); Caldweh 
et al. (2011) noted that Bo 82 was particularly 
red, hence the absorbing material is likely to reside 
in the cluster. The photon index T = 1.20±0.09; 
the unabsorbed 0.3-10 keV luminosity Lo.3-10 = 
2.6±0.2 X lO^s erg 8"^ and xVdof = 127/147. The 
observed photon index is lower than the expected 
range for black holes (F ~1.4-1.7); however, fixing 
r=1.4 resulted in a good fit with x^/dof <1; we 
present a x^/dof = 1 (x^ = 147) contour plot for 
r vs Nn in Fig. IH 

The best fit COMPTT model described a hot, op- 
tically thin corona, with kTo ^^54 keV and r ~1, 
for x^/dof = 127/145; however, these parameters 
could not be constrained, as kT^ lies well outside 
the XMM-Newton pass band. Indeed, electron 
temperatures of 100 or 300 keV (typical of low 
state black holes, see e.g. lMcClintock fc Remillard 
20061 and references within) were acceptable. 
Such spectra contrast with the cool, optically 
thick coronae obtained when modeling high lu- 
minosity neutron sta r binaries with COMPTT (e.g. 
Pi Salvo etallbOOOl ). 

The disk blackbody fit was acceptable (x^/dof 
= 140/147); however the inner disk temperature 
(kTin) was 3.2±0.4 keV, which is too hot for a 
black hole high state. Hence, we reject this fit on 
physical grounds, even though it is a statistically 
acceptable fit to the data. 

The two component model yielded A'h = 
4x10^^ and F = 1.2, as with the simple power 
law model. Furthermore, XSPEC was unable to 
estimate the uncertainties in the blackbody pa- 
rameters. We therefore conclude that any thermal 
component is too marginal to be detected; hence, 
XB082 is unlikely to be in the high accretion rate 
state for neutron star LMXBs. We present the 
spectrum and best fit two component model in 
Fig.lSl 

We conclude that XB082 exhibited non-thermal 
emission approximated by a power law. However, 
active galactic nuclei (AGN) also exhibit similar 
emission spectra, hence we calculated the proba- 
bility of a coincident AGN fr om the 2-10 keV lum i- 
nosity function provided bv lMoretti et al. (2003). 
The observed 2-10 keV flux of XB082 derived 
from the best fit power law model was 2.7x10"^^ 



bright existing within 1" of one of the 428 GCs 
is ~2xl0-^. We therefore conclude that XB082 is 
a candidate black hole LMXB. 

If we assume F=1.4, then the black hole would 
require a mass ~30 Mq for a 15-50 keV lumi- 
nosity <0 .1 J^Rnn. consiste nt with the low state 
(following iTang et al. 2010[ ). This is rather more 
massive than the Galactic stellar mass black holes 
( Ozel et allbOlOt ). but consistent with the primary 
in the dynamically confir med black hole -|- Wolf- 
Raye t binary ICIO X-1 ([Silverman fc Filippenko 
20081 ). 

We note that the observed spectrum of XB082 
(r ^1.2) is more often associated with high 
mass X- ray b i naries (HMXBs) with Be donors 
fsee e.g. iReie 201ll f or a recent review); e.g. 
Haberl fc PietschI (|20o3 ) examined XMM-Newton 
observations of 1 1 Be XBs in the Small Magellanic 
Cloud, and found the photon index distribution 
to be strongly peaked at F = 1.00, with a stan- 
dard devition of 0.16. Therefore XB082 could 
be a HMXB superposed on the globular cluster. 
However, we note that known Be XBs tend to be 
transient X-ray sources, as they have long, eccen- 
tric orbits and only accrete for a short time when 
the neutron star is closest to the donor; further- 
more, this emission spectrum is ge nerally con fined 



to the range -lO'^^'-lO-^^ erg s"^ ()Reigll201lL and 



references within). Persistently bright Be XBs 
exist, but tend to have X-ray luminosities <10'^^ 
erg s"^ If XB082 is a neutron star HMXB, it is 
unlike any know thus far. 

3.2. XB153 

XB153 was observed in 75 ACIS and 45 HRC 
observations, and appears to be persistently 
bright; its 0.3-10 keV luminosity varied by a factor 
^3 up to '^2.4x10'^^ erg s^^. It was also observed 
in the 2002 XMM-Newton observation with 60 ks 
good time. 

The best fit power law to the XMM-Newton pn 
spectrum yielded iVn = 8.5±1.0xl02° atom cm^^, 
and F = 1.62±0.04; xVdof = 429/444 and Lq. 



1.01±0.04xl038 erg s'^. 



erg cm 



the probability of an AGN this 



3-10 

Fitting a COMPTT 
emission model yielded a hot, optically thin corona 
with unconstrained kTo ~60 keV and t ^1. 

The disk blackbody model was rejected, with 
xVdof = 1067/445. 

Again, A^h and F for the two component model 
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were consistent with the values for the single power 
law, and XSPEC was unable to estimate the un- 
certainties for the blackbody parameters. We 
present the spectrum and two component fit in 
FigH 

XB153 also appears to have been in a non- 
thermal state. The observed 2-10 keV flux for the 
best fit power law model was 8.6x 10"^'^ erg cm ~^ 
s~^; hence, the probability of finding an AGN of 
this brightness within 1" of any of the 428 GCs is 
1.3x10"'^. The luminosity for the XMM-Newton 
observation was a factor ~2 lower than the peak. 

3.3. XB163 

XB 163 was observed in 10 ACIS and 29 HRC 
observations. XB163 exhibited at least 5 out- 
bursts over the ~4000 day viewing period, and 
was brightest during ACIS observation 8184 (2007 
February 14); the exposure was 5 ks. The 0.3- 
10 keV luminosity of XB163 varied by a factor 
^^650, hence it cannot be an AGN. XB163 was 
observed three times in the public XMM-Newton 
archive; once in 2006, December and twice in 2008, 
February; however, one of the 2008 observations 
experienced severe flaring. We have added the 
two flare-free XMM-Newton observations to the 
XB163 lightcurve in Fig. Et. The first XMM- 
Newton observation was 45 days before the bright- 
est Chandra observation, and may be part of the 
same outburst. The second XMM-Newton obser- 
vation was during an outburst that was missed 
by Chandra; although an ACIS observation was 
conducted 37 days beforehand, the roll angle was 
unfavorable. 

The best spectrum was obtained from Chandra 
observation 8184, as XB163 did not significantly 
exceed the neutron star luminosity threshold for 
the low state in the XMM-Newton observations. 
The mean intensity of XB163 during this obser- 
vation was ~0.1 count s~^, which would result 
in significant pileup for an on-axis X-ray source. 
However, XB163 was ~14' off- axis during this ob- 
servation. Each photon is expected to be detected 
in a 3x3 array of ACIS pixels; the brightest 3x3 
portion of the XB163 image yielded only 17 pho- 
tons, equivalent to one photon per ~100 frames; 
hence, the pileup was negligible. 

The best fit power law model, with photon in- 
dex F = 1.5±0.2, and TVh = l.ltar xlO^° atom 



cm ^, provided a good fit, with x^/dof =26/22; 
the 0.3-10 keV luminosity for this best fit power 
law model was 1.40±0.13x lO^^ erg s'^ The 
COMPTT model again yielded a hot, thin corona, 
with unconstrained parameters: kTo ^52 keV, r 
-1. 

The disk blackbody fits required A^h to be 
frozen at 6.7x10^° atom cni^^, to prevent the 
absorption from falling below the Galactic line 
of sight absorption; kT = 1.4±0.2 keV, con- 
sistent with the temperature range exhibited 
by black hole binaries i n thei r high state (see 
iMcClintock fc RemiUardl 120061 and references 
within). The 0.3-10 keV luminosity for this fit 
was 9.8±1.8xl03^ erg s'^ 

Initial fitting of the two component model fa- 
vored a temperature of 199 keV and normalization 
~0, with all the observable flux coming from the 
power law, for x^/dof =27/21. However, a better 
fit was found with a 0.7±0.3 keV blackbody and 
a F = 1.6±0.4 power law, giving x^/dof = 24/21; 
the absorption had to be frozen at the Galactic 
line of sight column density. The spectrum and 
best two component fit is presented in Fig[3l 

We cannot determine whether the accretor is a 
neutron star or black hole from the spectra alone. 
However, the highest observed 0.3-10 keV lumi- 
nosity is 1.4x 10"^^ erg s~^, a factor ~2 higher than 
the peak brightnesses of any of the >20 outbursts 
observed from the Galactic neutron star tra nsients 
Aql X-1 and 4U 1608-52 (|Lin et all 120071) . Fur- 
thermore, XB163 was active ~45 days before our 
observed peak, with a 0.3-10 keV luminosity of 
2.3±0.6xl03^ erg s"\ with 90% confidence un- 
certainties; if this activity is related to the out- 
burst, then the true peak may have been con- 
siderably brighter. We note that the M31 tran- 
sient CX0M31 J004253. 1-1-411422, discovered in 
2009 with a 0.2-10 keV luminosity -4x10^^ erg 



( Henze et al.ll200^ . maintained a luminosity 



^lO"^** erg s~ 150 days after discovery (Nooraee 
et al, in prep); hence, the true peak luminosity 
may have been considerably higher than the ob- 
served peak. We find it unlikely that XB163 is a 
high accretion rate neutron star binary, and sug- 
gest that it is a black hole binary. However, if the 
outburst detected in the XMM-Newton observa- 
tion were unrelated to the Chandra peak, then a 
neutron star primary would still be feasible. 
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3.4. XB185 

XB185 was observed in 25 ACTS and 45 HRC 
observations due to its high off-axis angle, but ap- 
pears to be persistently bright. The 0.3-10 keV 
luminosity varies by a factor ~3, up to ^^IxlO'^^ 
erg s~^. It was also observed in the 60 ks XMM- 
Newton observation. 

The XMM-Newton pn spectrum was well 
described by a power law model with Nn = 
1.09±0. 17x10^1 atom cm'^, and T = 1.64±0.05; 
X^/dof = 242/232. Modeling the emission with a 
COMPTT model yields kT^ ~60 keV and t ^0.9; 
again these parameters are unconstrained. 

The disk blackbody fit was rejected with x^/dof 
= 405/233. Once more, A'h and F are consistent 
with the power law model for the two component 
fit, while the blackbody parameters could not be 
constrained; the spectrum and best two compo- 
nent fit is presented in Fig. [31 

The 2-10 keV observed fiux for the best fit 
power law model was 5.9x10"^'^ erg cm^'^ s~^. 
The probability for an AGN of equivalent bright- 
ness being within 1" of one of the 428 GCs is 
2x10^5. 

3.5. Properties of the host clusters 

We present the I magnitude, FK colours, 
masses ages and [Fe/H] metalicities of the four 
GCs suspected of hosting new BHCs in Table [2l 
along with the two M31 GC BHCs already identi- 
fied; the magnitudes were obtained from the RBC; 
the masses, ages and metal icities were drawn from 
Caldwell et all (l2009h and ICaldwell et al.l (l201lh . 
Comparison with the prop erties of the general GC 
population presented by IPeacock et aD (|2010at ) 
suggests that the GCs harboring BHCs are more 
massive (brighter), and redder / more metal rich 
than the general population. We present in Fig. [S] 
the I vs. FK color magnitude diagram for M31 
GCs; points represent the whole GC population 
with I and K magnitudes, open circles represent 
GCs in our field with associated X-ray emission 
(see Barnard et al., in prep), and filled circles 
represent the six M31 GCs harboring BHCs; the 
mean of the whole GC population is represented 
by a star, the mean of the non-BHC X-ray GCs 
by a triangle, and the mean of the BHC GCs is 
represented by a square. It appears that the BHC 
GCs are a rather massive and red subset of the 



GCs associated with X-ray sources. 

While all of the BHC GCs are 2-14 times more 
metal rich than the mean M31 GC m etalicity of 
-1.08 found bv ICaldwell et al.l (|201l[ ). and more 
massive than 63% of the GC population, only Bo 
163 is particularly massive and metal rich. Bo 82 
is the 4th most massive of th e 379 GCs analysed 
bv lCaldwell et al.l(|2009ll201l[ ). but only the 102nd 
most metal rich. Bo 144 is the 139th most mas- 
sive, but the 10th richest; indeed, i t is supersolar 



and r icher than any Galactic GC ([Caldwell et al 
201ll ). Bo 153 is the 87th most massive, but the 
36th richest. Bo 163 is the 20th most massive and 
the 21st richest. Bo 185 is the 53rd most massive 
and the 75th richest. Hence GCs that are either 
massive or metal rich are able to produce bright 
X-ray sources, as well as GCs that are both. 

4. Discussion 

Despite early indications for an absence of stel- 
lar mass black hole binaries in globular clusters, 
they are becoming increasingly common. Out of 
the 35 X-ray sources associated with globular clus- 
ters in the central region of M31, 5 harbor BHCs. 
Four of these appear to be persistently bright, 
and ar e consistent with the theoretical predic- 
tions of Kalogera et al. ( 20041) for binaries formed 
by tidal capture of a main sequence star, or for 
ultra-compact black hole -I- white dwarf binaries 
( Ivanova et al. 2010l ). However, XB163 is a recur- 
ring transient, and provid es the first test of the 
theoretical predictions of iKalogera et al. (2004) 
regarding binaries formed by exchange. We find 
that the GCs that are metal rich or massive are 
able to produce bright X-ray sources, in addition 
to GCs that are both. 

We observed 5 outbursts in the Chandra and 
XMM-Newton observations of XB163 over -4000 
days; however, the large off-axis angle for XB163 
in these observations means that further out- 
bursts may have been missed. Four of the out- 
bursts occurred over 1100 days. Furthermore, 
Trudolvubov fc Priedhorskvl (|2004l ) found 3 out- 
bursts within 400 days in the ROSAT observa- 
tions. Such behavior may be consistent with a 
black hole -I- main sequence star formed by ex- 
change |Kalogeraet_alJ (|2004[ ). or due to the com- 
plex behaviour of a black hole -I- white dwarf bi- 
nary in a hierarchical triple system, as envisioned 
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bv llvanova et al. I (|2010l) . 

We shall now discuss alternative explanations 
for the observed behaviour. We shall discuss 
whether they are coincident AGNs, blended emis- 
sion from multiple sources, or neutron star binaries 
with beamed emission. 

As discussed in the previous section, the prob- 
ability of each source being an AGN that is coin- 
cident with one of the GGs in our field is already 
small. The probability that XB082, XB153, X185 
and XB144 are all coincident AGNs within 1" of 
a GC is <10~^'~', since the X-ray luminosities ob- 
served in the XMM-Newton observations used to 
obtain the spectra were each lower than the peak 
observed luminosity in the Chandra data. 

We shall now consider the possibility that the 
BHCs in our field are instead blends of multiple 
X-ray sources; 4-14 1.4 Mq neutron star bina- 
ries accreting at 10% Eddington would be required 
to produce the observed XMM-Newton spectra 
of XB082, XB153 and XB185, and most would 
have to be persistent to account for the Chandra 
hghtcurves. We find that XB082 varied by 10^* 
erg in 9 days, XB153 varied by 10^^ erg s~^ in 
2 days, and XB185 varied by 3 x 10^'' erg s^^ 
over ^ 2 hours. Such variation is more likely to 
come from a single source than the concerted vari- 
ation of the several low state neutron star binaries 
require to matched the observed luminosities. 

Finally we discuss the possibility that these 
BHCs are neutron star X-ray binaries with 
beamed emission. To do this, we estimate the 
minimum beaming factor, b, which we define as 
the unabsorbed 15-50 keV luminosity for each 
source, divided by 0.1 LEdd for a 2 Mq neutron 
star (2.6x10'^^ erg s~^); since X-ray binaries in 
the low state have electron temperatures '^100- 
300 keV, we may assume power law emission when 
converting to 15-50 keV flux. The probability that 
we intercept beamed radiation from a randomly 
oriented X-ray binary system is then < 1/b. As- 
suming the 90% confidence lower limits to A^h and 
r (i.e. obtaining the lower limit to the 15-50 keV 
luminosity) yields a probability < 4 x 10^^ that 
these are neutron star LMXBs with beamed emis- 
sion. If we include XB045 (Barnard et al. 2008), 
then the probability that they are all beamed neu- 
tron star LMXBs falls to <3xl0-5. 

In total, we have identified 5 BHCs in M31 GCs 



that a.ppear to be persistently bright (jBarnard et al 
20081 : iBarnard fc Kolbll2009( . and this work); 4 of 
these are in the central region of M31, and the 
null hypothesis that AGNs coincide within 1" of 
4 out of 428 GCs in our field has a probability of 
<lxl0^^°, using the 2-10 kcV luminosity func- 
tion for AGNs provided by iMoretti et all (i2003h 



Before these recent discoveries, iKalogera et al 
(|2004l ) inferred that tidal capture of main sequence 
stars by black holes resulted in the disruption of 
the main sequence star. It is now increasingly 
likely that either tidal capture X-ray binaries can 
survive, or our BHCs could b e ultracompact Blac k 
hole -I- white dwarf binaries ( Ivanova et al. 2010f l. 
We therefore suspect that other GC BHCs have 
been already observed, but not yet recognized. 

We thank the referee for their insightful com- 
ments; this work was significantly improved as 
a result. This research has made use of data 
obtained from the Chandra data archive, and 
software provided by the Chandra X-ray Center 
(CXC). This work also used observations made 
by XMM-Newton, an ESA science mission with 
instruments and contributions directly funded by 
ESA member states and the US (NASA). R.B. is 
funded by Chandra grant GO9-0100X and HST 
grant GO-11013. M.R.G. is partially supported 
by NASA grant NAS-03060. We thank Andrew 
King for useful discussions. 
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Fig. 1. Merged ACIS image of the central region of M31, with radius ~20', in the 0.3-7.0 keV band; our 
BHCs are circled and labeled. 
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Fig. 2. — Long-term 0.3-10 keV luminosity lightcurves from up to 78 ACIS and 45 HRC observations, 
indicated by black crosses and red circles respectively; blue diamonds represent XMM-Newton observations 
of XB163. For each lightcurve we provide the best fit line of constant intensity Leon- The four lightcurves 
are: (a) XB082— Lco„ = 2.33±0.05xl038 erg s'^, x^/dof = 181/47; (b) XB153— Leon =1.029±0.009xl038 
erg s-\ x^/dof = 817/119; (c) XB163— Leon = 1.5±0. 2x10^6 erg s^S x^/dof =680/40; (d) XB185~ Leon 
= 5.61±0.09xl0'^^ erg s^^, x^ /dof — 185/69. A color version is available in the electronic edition. 
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Fig. 3. — Unfolded spectra multiplied by energy, plotted with the best fit two component model, consist- 
ing of a blackbody (red dots) and a power law (black dots): (a) XB082, from XMM- Newton observation 
0109270101 with 3110 net source counts; (b) XB153 from XMM-Newton observation 0112570101 with 11923 
net source counts; (c) XB163 from ACIS observation 8184 with 530 net counts; (d) XB185 from XMM- 
Newton observation 0112570101, with 5345 net source counts. The blackbody components in XB082, XB153 
and XB185 were too marginal to be constrained. A color version is available in the electronic edition. 
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Fig. 5. — Detail ol the I vs I-K band colour diagram for the M31 GCs with available data in the REG; 
some outliers are excluded for clarity. The general population is represented by dots, GCs with associated 
X-ray emission with open circles, and GCs harbouring our BHGs with filled circles. The mean colours and 
magnitudes of the total GC population. X-ray GCs with BHGs and X-ray GCs without BHCs are represented 
by a star, a square and a triangle respectively. A color version is available in the electronic edition. 
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Table 1: Spectral properties of the four black hole candidates. For each source, we give the observation used 
to get the best spectrum, and the net source counts. We then give the absorption / 10^^ atom cm^, the 
photon index and x^/dof for the best fit power law model. Numbers in parentheses indicate 90% confidence 
uncertainties on the last digit. 



Source 


Obs 


Counts 




r 


xVdof 


7-37 


XB082 


XMM 0109270101 


3110 


3.9(3) 


1.20(9) 


127/147 


26(2) 


XB153 


XMM 0112570101 


11923 


0.85(10) 


1.62(4) 


429/444 


10.1(4) 


XB163 


Chandra 8184 


567 


1.0{tV) 


1.5(2) 


26/22 


14(3) 


XB185 


XMM 0112570101 


5345 


1.09(7) 


1.64(5) 


242/232 


7.0(5) 
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Table 2: Spectral properties of the six M31 glo bular clusters suspected of harboring b l ack h ole candidates, 
including two that were previously identified (jBarnard et al. 2008t Barnard fc Kolb 2009( ). Magnitudes 
were obtained f rom t he R evised Bologna Catal ogue, while the age, mass and metalicities were taken from 
Caldwell et all (|20n9l) andlCaldweh et alJ (|201lh. 



GC 


1 


I-K 


Log Age / yr 


Log Mass / M0 


[Fe/H] 


Bo 45 


14.51 


1.60 


10.2 


5.94 


-0.86 


Bo 82 


13.63 


2.18 


10.0 


6.72 


-0.74 


Bo 144 


15.29 


2.23 


10.1 


5.58 


+0.08 


Bo 153 


14.94 


1.84 


10.1 


5.81 


-0.28 


Bo 163 


13.71 


2.02 


10.1 


6.26 


-0.13 


Bo 185 


14.36 


1.72 


10.1 


6.03 


-0.61 
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